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Trigonometric Parallaxes of Massive Star Forming Regions: III. 

G59.7+0.1 and W 51 IRS2 

Y. Xu^'^ M. J. Reid^ K. M. Menten\ A. Brunthaler\ X. W. Zheng^, L. Moscadelli^ 

ABSTRACT 

We report trigonometric parallaxes for G59. 7+0.1 and W 51 IRS2, corre- 
sponding to distances of 2.161q;o9 ^P"^ ^-^-ia ^pc, respectively. The distance 
to G59. 7+0.1 is smaller than its near kinematic distance and places it between 
the Carina-Sagittarius and Perseus spiral arms, probably in the Local (Orion) 
spur. The distance to W 51 IRS2, while subject to significant uncertainty, is 
close to its kinematic distance and places it near the tangent point of the Carina- 
Sagittarius arm. It also agrees well with a recent estimate based on 0-type star 
spectro/photometry. Combining the distances and proper motions with observed 
radial velocities gives the full space motions of the star forming regions. We find 
modest deviations of 5 to 10 km s~^ from circular Galactic orbits for these sources, 
both counter to Galactic rotation and toward the Galactic center. 

Subject headings: techniques: interferometric — astrometry — galaxy: structure 
— stars: individual (G59.7+0.1, W 51 IRS2) 

1. Introduction 

We are carrying out a large project to study the spiral structure and kinematics of 
the Milky Way by measuring trigonometric parallaxes and proper motions of star forming 
regions. The target sources are 12 GHz methanol masers and we use the National Radio 
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Astronomy Observatory ' s B Ve ry Long Baseline Array (VLBA). Details of this prog: 
be found in iReid et al. I (120081 ) . hereafter called Paper I. 



ram can 



Here we present observations of G59.7+0.1 (IRAS 19410+2336) and W 51 IRS2. De- 
pending on its distance, G59. 7+0.1 could be in either the Carina-Sagittarius or Local spiral 
arm of the Milky Way. W 51 IRS2 is a very well studied region of high-mass star formation; 
its radial velocity is very close to the maximum allowed by circular rotation for standard 
models of the Milky Way, and it has generally been assigned the tangent-point distance, 
D, in the Carina-Sagittarius arm. Our direct measurements of distance and proper motion 
indicate the true location of these star forming regions in the Galaxy and their departures 
from circular Galactic orbits. 



Observations and Calibration 



We conducted phase-referenced observations of G59. 7+0.1 and W 51 IRS2 with re- 
spect to two extragalactic radio sources with the VLBA under program BRIOOD in order to 
measure parallaxes. Paper I describes the general observational setup and data calibration 
procedures, so here we only describe details specific to the observations of G59. 7+0.1 and 
W 51 IRS2. 

The time between epochs was planned to be three months, matching the eastward and 
northward extrema of the Earth's orbit as seen by the sources. The observations were 
conducted on 2005 Jul. 13 and Oct. 20, 2006 Jan 15, Apr. 23 and Oct. 19, and 2007 Apr. 
19. However, the data for the epoch of 2006 Apr. 23 were not correlated at the position of 
W 51 IRS2 and were lost for this source. 

Background compact extragalactic sources were chosen as follows. For G59. 7+0.1, we 
selec ted J1946+2300 (with separation of 1.0° from the maser target) from the ICRF source 
list jMa et al. \\l99^ and J1941+2307 (separation 7°) and J1943+ 2330 (separation 0.3°), 
based on a VLA survey of compact NVSS sources (IXu et al. Il2006al ). Ultimately, we used 
only data from the first two sources, as we failed to detect J1943-I-2330 . For W 51 IRS2, we 
chose two sources from the VLBA Calibrator Survey (jPetrov et al. II2006I ). J1922+1530 (sepa- 
ration 1.0°) and J1924+1540 (separation 1.2°), augmented by J1922+1504 (separation 0.6°), 
which was found in our VLA survey. Two strong sources (J1800+3848 and J1922+1530) 
were observed near the beginning, middle and end of the observations in order to monitor 
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delay and electronic phase differences among the IF bands. 

After applying the basic calibration procedures described in Paper I, we used the maser 
features toward G59. 7+0.1 at Vlsr = 27.4 km s~^ and toward W 51 IRS2 at Vlsr = 
56.4 km s~^ for interferometer phase-reference data. When imaging the data referenced 
to G59. 7+0.1, we adopted a round restoring beam of 1.7 mas (FWHM), slightly larger than 
the geometric mean of the interferometer response ("dirty beam") of 1.7 x 1.1 mas at a posi- 
tion angle of 14° East of North. For the data referenced to W 51 IRS2, we adopted a round 
restoring beam of 3.3 mas (FWHM), slightly larger than the dirty beam of 3.3 x 1.5 mas at 
a position angle of 134°. 



3. Parallax and Proper Motion 
3.1. G59.7+0.1 

A map of the maser emission made by integrating all channels with detectable emission 
at the first epoch is shown in Fig. [TJ Millimeter and centimeter wavelength continuum data 
associated with the maser are also presented. We analyzed 8.4 GHz data, which we obtained 
from the VLA archival database, and find a weak continuum source with flux density of 
0.77 + 0.18 mJy at (aj2ooo, (^J2ooo) = (19 43 11.21, +23 44 03.32), with a probable position 
uncertainty of about 0.2 arcsec. We conclude that the methanol maser (Table 1) and the 
compact continuum sources are associated. The continuum source has an upper limit of 
^ 2" for its size. One milli meter wavelength core is also asso ciated with the maser and the 



weak cm continuum source (IBeuther. Schilke fc Stanke 1120031 ). 



We show the first epoch images of each background continuum source in Fig. [2l One 
can see that they both appear dominated by a single compact component. 

When conducting phase-referenced observations, it is important that the position of 
the reference source matches the interferometer phase center in order to minimize second- 
order positional errors and improve image quality. Since the ICRF source J1946+2300 has 
a position accuracy of ~ 1 mas, we used its position as the basis for all absolute positions 
given in Table [H 

We fitted elhptical Gaussian brightness distributions to two maser spots and the two 
background radio sources for all seven epochs. In Fig. [3l we plot the positions of two 
maser spots (at Vlsr of 26.6 and 27.4 km s~^) relative to two background sources. The 
measured positions of the G59. 7+0.1 masers were then modeled as a linear combination 
of the elliptical parallax and linear proper motion signatures. Because systematic errors 
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Fig. 1. — Upper Panel: Continu um emission at 250 GHz (contour s) imaged with the IRAM 
Plateau de Bure interferometer (IBeuther. Schilke fc Stanke II2003I ). The the star marks the 
12.2 CH3OH maser position (see Table 1). Middle Panel: Weak 8.4 GHz emission toward 
G59. 7+0.1 imaged with the VLA. Contours represent 30, 50, 70, and 90% of the peak 
brightness of 0.68 mJy beam^^. The lowest contour represents 2 times the rms noise level, 
while the thick 50% contour closely follows the size of the synthesized beam (2.9" x 2.4"). 
The (0,0) position corresponds to the position of the 12 GHz methanol maser. Bottom 
panel: Map made by integrating all channels with detectable methanol maser emission at 
the first epoch. The li,sR = 27.4 km s~^ maser spot, located at the origin, was used for 
the parallax fits. The maser spot located near (—0.001, —0.005) is at Flsr = 26.6 km s~^. 
Contours represent 30, 50 (bold), 70, and 90% of the velocity-integrated flux density of 0.82 
Jy km s~^. The restoring beam is 1.7 mas FWHM for the maser image. 
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Fig. 2. — First epoch images of the two background sources phase-referenced to the 
G59. 7+0.1 methanol maser. For both sources the lowest negative (dashed) and positive 
(solid) contours are 4 times the rms noise level, which is 0.4 mJy beam~^ for J1941+2307 
(upper panel) and 0.7 mJy beam~^ for J1946+2300 (lower panel). Other contours repre- 
sent increases by a factor of 2 starting from that level. The peak brightness and integrated 
flux density arc 10.8 mJy beam^^ and 16.3 mJy, respectively, for J1941+2307 and 74.8 
mJy bcam^^ and 114 mJy, respectively, for J1946+2300. The upper limit of source sizes is 
3.0 mas for J1941+2307 and 2.9 mas for J1946+2300. The FWHM of the synthesized beams 
are represented in the lower left corner of each panel. 
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(owing to small uncompensated atmospheric delays and, in some cases, varying maser source 
structure) typically dominate over signal to noise considerations when measuring relative 
source positions, we added "error floors" in quadrature to the formal position uncertainties. 
We used different error floors for the Right Ascension and Declination data and adjusted 
them to yield post-fit residuals with per degree of freedom near unity for both coordinates. 
Individual fits are given in Table 2. 

Fitting for the parallax and proper motion for both sources simultaneously, we obtain 
71 = 0.463 ± 0.020 mas. The quoted parallax uncertainty is the formal fitting uncertainty, 
multiplied by to account for possible correlations between the position data for the two 
maser spots. This parallax corresponds to a distance of 2.16lg o9 ^P^' which is smaller than 
the "near" kinematic distance of 2.7 kpc, and rules out the far distance of 5.8 kpc. The 
average proper motions in the eastward and northward directions are —1.65 ± 0.03 and 
—5.12 ± 0.08 mas y~^, respectively, as listed in Table [2l Similarly, the uncertainties were 
also multiplied by -\/2. 



3.2. W 51 IRS2 

We made a 12 GHz methanol maser map of W 51 IRS2 by integrating all channels 
with detectable emission at the first epoch. This map is shown in Fig. HI superposed on a 
continuum (23.8 GHz) image from archival VLA data (AS724). 

In Fig. [5], we show maps of the background continuum sources, phase referenced to the 
Vlsr = 56.4 km s~^ maser spot in W 51 IRS2, from the first epoch. These background 
sources are dominated by a single compact component. Absolute positions for the maser 
reference spot and the background sources are given in Table [H These ar e based on the 



position of J1924+1540, which is uncertain by ~ 1 mas (jPetrov et al. II2006I ) 



In order to determine the parallax and proper motion of W 51 IRS2, we used positions 
of two strong maser spots relative to all three background sources. Following the fitting pro- 
cedures discussed for G59. 7+0.1, individual spot/background source parallax solutions are 
listed in TableO A combined parallax solution yielded vr = 0.195±0.071 mas, corresponding 
to a distance of 5.l!l^^4 kpc, which will be discussed in §3.31 The data and model used for this 
fit are shown in Fig. El The formal parallax uncertainty has been multiplied by \/2, in order 
to account for possible correlations between the position data for the two maser spots. Some 
of the relative positions for the fifth epoch (2006.80) appear to be outliers. Were we to drop 
the data from this epoch, the parallax estimate decreases to tt = 0.166 ±0.069 mas, suggest- 
ing a slightly greater distance. The average proper motions of the two maser spots in the 
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East Offset (mas) Epocll (years) Epoctl (years) 

Fig. 3. — Parallax and proper motion data and fits for G59. 7+0.1. Plotted are position 
measurements of two maser spots at Vlsr = 27.4 and 26.6 km (open and solid symbols) in 
G59. 7+0.1 relative to two background sources: J1941+2307 (red triangles) and J1946+2300 
(green squares). Left Panel: Positions on the sky with first and last epochs labeled. The 
expected positions from the parallax and proper motion fit are indicated (crosses). Middle 
Panel: Eastward (solid lines) and northward (dashed lines) positions and best fit parallax 
and proper motions fits versus time. Right Panel: Same as the middle panel, except the best 
fit proper motions have been removed, allowing all data to be overlaid and the effects of only 
the parallax seen. 
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Fig. 4. — The large image shows weak 23.8 GHz continuum emission from a hypercompact 
HII region south of W 51 d, the ultracompact HIT region associated with W 51 IRS2. The 
image was produced from archival VLA with a restoring beam (indicated in the lower left 
corner) of ~ 0''08. Contours are -5, 5, 10, and 20 times the rms noise of 0.28 mJy beam~^. 
The methanol maser is, as indicated, clearly associated with the hypercompact HII region. 
The separation between them is around 0.07 arcsec. The inset shows the velocity-integrated 
methanol maser emission with contours representing 5, 10, 20, 40, 80, and 160 times 5 
mJy beam~^ km s~^. The x- and y-axes give east and west offset, respectively, relative to 
the position given in Table 1. 
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Fig. 5. — First epoch images of of the three background sources phase-referenced to the W 
51 IRS 2 methanol maser. For all sources the lowest negative and positve contours are 4 times 
the rms noise level, which is 0.6 mJy beam"^ for J1922+1530 (top panel), 1.8 mJy beam~^ 
for J1922-M504 (middle panel) and 2.9 mJy beam ^ for J1924-I-1540 (bottom panel). Other 
contours represent increases by a factor of 2 starting from that level. The peak brightness and 
integrated flux density are 22.5 mJy bearn"^ and 43 mJy, respectively, for J1922-I-1530, 24 
mJy beam~^ and 37 mJy, respectively, for J1922-I-1504, and 465 mJy beam~^ and 615 mJy, 
respectively, for J1924-I-1504. The upper hmit of source sizes is 2.6 mas for J1922-I-1530, 3.8 
mas for J1922+1504, and 0.0 mas for J1924+1504. The FWHM of the synthesized beams 
are represented in the lower left corner of each panel. 
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eastward and northward directions are —2.49 ±0.08 and —5.51 ±0.11 mas y ^, respectively, 
as listed in Table [31 



3.3. The Distance to W 51 IRS 2 



It is interesting to compare our directly measured distance of 5.1^ ?'^ kpc with other 



distance determinations that have recently published for W51 IRS 2. Ilmai et al.l (|2002| ) 
performed multi-epoch VLBI observations of the intense H2O maser outflow in the region 
("W 51 N"). A kinematic model they fitted to the measured internal proper motions of the 
outflow contains its distance as free parameter, for which they find a best fit value of 6.1±1.3 
kpc. 



Barbosa et al.l (120081 ) combine near infrared spectroscopy and radio continuum and re- 



combination line observations to characterize the exciting star of W 51 IRS 2. Using available 
data and calculations modeling the temperatures, radii and Lyman continuum production 
rates of Zero Age Main Sequence (ZAMS) stars, they propose a classification of its spectral 
type as 03 or 04 (ZAMS). Using the source's observed bolometric luminosity, they derive 
D = 5.1 and 5.8 kpc for these two choices, respectively, which compares very well with our 
distance. 

We note that the same group recently also obtained a distance of 2.0 ± 0.3 from spec- 
troscopic and photomet ric observations of four objects in W 51 A classified as 0-type stars 
(IFigueredo et al. I l2008l ). W 51 A and W 51 IRS 2, separated by an angle of roughly 1 
arc minute (or 1.5 pc), are commonly thought to belong to the same complex and a (near) 
kinematic distance of 5.5 kpc has been found by iRusseill (120031 ) for W 51 A, comparable to 
our parallax distance of IRS 2. 

Currently, large-scale infrared surveys are leadi ng to the discovery of many new open 



star clusters thr oughout the Milky Way (see, e.g., Ilvanov et al.l l2002l : jPutra et al.l 12003 
Bica et al.l l2003l ). why may trigger a renaissance in efforts to use spectro/photometry of 
0-ty pe stars for distance estimates with the goal of constraining Galactic structure (see, 
e.g., iMessineo et al.ll2006l ). In view of this, we find that the discrepancy discussed above 



deserves further investigation. We note that in the famous case of W30H in the Perseus arm 
the -star method yielded D = 2.3 kpc for the close-by Per OB 1 association ([Humphreys 
19781 ). roughly half the kinematic distance implied by a model of Galactic rotation. This 
turned out to be in excellent agreement with the 2.0 kpc directly determin ed via two VLBI 



paral l ax measurements of CH3OH and II2O masers in W30H much later ([Hachisuka et al. 



20061 ; IXu et al. Il2006bl ). Future maser VLBI trigonometric parallax measurements of W 51 
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Fig. 6. — Parallax and proper motion data and fits for W 51 IRS2. Plotted are position 
measurements of two maser spots at Vlsr = 56.0 and 56.4 km (open and solid symbols) 
relative to the three background sources: J1922+1504 (blue hexagons), J1922+1530 (red 
triangles) and J1924+1540 (green squares). Left Panel: Positions on the sky with first and 
last epochs labeled. Data for the two maser spots are offset horizontally for clarity. The 
expected positions from the parallax and proper motion fit are indicated (crosses). Middle 
Panel: Eastward (solid lines) and North (dashed lines) positions and best fit parallax and 
proper motions fits versus time. Data for the two maser spots are offset vertically and small 
time shifts have been added to the data for clarity Right Panel: Same as the middle panel, 
except the with the best fit proper motions have been removed, allowing all data to be 
overlaid and the effects of only the parallax seen. 
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A could certainly contribute to the perplexing question described above. 



3.4. Galactic Locations and 3-D Motions 



In order to study the 3- dimensional motion of the maser sources in the Galaxy, we 
converted the proper motions and radial velocities to a Galactocentric reference frame. We 
adopt the lAU standard constants of Rq = 8.5 kpc and Gq = 220 km s~^, and the Hipparcos 
Solar Motion values U = 10.0 ± 0.40, V = 5.25 ± 0.60, and W = 7.17 ± 0.40 km s"^ from 



Dehnen fc Binney I (119981 ). For these values and a flat rotation curve for the Galaxy, the 
Galactocentric distance of G59. 7+0.1 is 7.64 kpc. Its velocity in the direction of Galactic 
rotation is 10 ± 3 km s~^ slower than for a circular orbit. Its velocity toward the Galactic 
Center is 7 ± 1 km s~^, and its velocity toward North Galactic Pole is —4 ± 1 km s~^. 
These uncertainties include measurement errors, but do not include systematic terms from 
uncertainty in Rq and Gq. Thus we find that G59. 7+0.1 has a peculiar motion of ~ 12 km s~^ 
directed mostly counter to Galactic rotation and toward the Galactic Center. 

Our trigonometric parallax places G59. 7+0.1 in the Milky Way between the Carina- 
Sagittarius and Perseus spiral arms. It seems to be a fairly distant member of the Local 
(Orion) arm or spur, located close to the point where the spur joins the Carina-Sagittarius 
arm. G59. 7 +0.1 is near the open cluster NGC 6823, wh ich is also thought to be located in the 



Local arm (Basharina. Pavlovskava fc Fill 



serve d for many galaxies (lAalto et al. Ill999[ Scoville et al. II2001I : iLa Vigne. Vogel fc Ostriker 
20061 ) . These spurs may form as a consequence of gravitational instabilities inside spiral arms 



pova 



19801). Spur-hke structures have been ob- 



or/an d effects of magne t ic fie lds (IBalbus I Il998l : iKim fc Ostriker 1 120021 : IShetty fc Ostriker 



20061 ). iKim fc Ostriker I (120021 ) showed that the growth of spurs can occur due to the mutual 
contributions of self-gravity and magnetic fields via the so-called magneto- Jeans instability. 

Adopting a distance of 5.1 kpc from the Sun, places W 51 IRS2 6.5 kpc from the Galactic 
center, in the the Carina-Sagittarius spiral arm and reasonably close the spiral arm tangent 
point. Converting the proper motions and radial velocity of W 51 IRS2 to a Galactocentric 
reference frame, we find a velocity in the direction of Galactic rotation that is 5 ± 10 km s~^ 
slower than for a circular orbit. Its velocity toward the Galactic Center is 21 + 15 km s~^, and 
its velocity toward North Galactic Pole is —3 + 5 km s^^. These peculiar motion uncertainties 
are fairly large, primarily because of the uncertain parallax measurement. 



4. Conclusions 



We have measured the parallax and proper motion of methanal masers in two regions of 
high-mass star formation. G59. 7+0.1 lies at a distance of 2.I6IQJ9 ^P*^ ^^'^ Local (Orion) 
arm or spur. Its space motion, relative to a frame rotating with the Milky Way, is about 
12 km counter to Galactic rotation and toward the Galactic center. Our parallax for 
W 51 IRS2, while less accurate than for G59.7+0.1, indicates that W 51 IRS2 is in the 
Garina-Sagittarius spiral arm. 

This work was supported by Ghinese NSF through grants NSF 10673024, NSF 10733030, 
NSF 10703010 and NSF 10621303, and NBRPG (973 Program) under grant 2007GB815403. 
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Table 1. Positions and Brightness 



Source R.A. (J2000) Dec. (J2000) (f) Brightness Vlsr Restoring Beam 

') (° ' ") (°) (Jy/beam) (km s"^) (mas, mas, deg) 

G59.7+0.1 19 43 11.2470 23 44 03.315 1.4 27.4 1.7 

J1946+2300 19 46 06.25140 23 00 04.4145 1.0 0.073 2.2x3.2 @ 33 

J1941+2307 19 41 55.1114 23 07 56.525 0.7 0.010 2.3x3.8 @ 30 

W51 19 23 39.8244 14 31 04.953 2.2 56.4 3.3 

J1924+1540 ... 19 24 39.45588 15 40 43.9417 1.2 0.46 2.4x4.3 @ -51 

J1922+1530 ... 19 22 34.6993 15 30 10.0327 1.0 0.23 1.9x4.2 @ -48 

J1922+1504 ... 19 22 33.2728 15 04 47.537 0.6 0.022 1.9x4.6 @ -43 



Note. — (p is the separations. The radial velocity of the masers and the size and shape of the 
interferometer restoring beam are listed for the first epoch's data. The position angle of the beam is 
defined as East of North. 

Table 2. G59. 7+0.1 Parallax & Proper Motion Fit 



Maser Vlsr Background Parallax fx^ 



(km s ^) Source (mas) (mas y ^) (mas y ^) 



26.6 J1941+2307 0.459 ±0.043 -1.68 ±0.06 -5.18 ±0.12 

26.6 J1946+2300 0.484 ± 0.025 -1.68 ±0.03 -5.09 ±0.10 

27.4 J1941+2307 0.436 ± 0.036 -1.63 ±0.04 -5.17 ±0.11 

27.4 J1946±2300 0.466 ± 0.027 -1.63 ±0.04 -5.08 ± 0.08 

26.6 combined 0.463 ± 0.020 -1.68 ± 0.03 -5.13 ±0.08 

27.4 -1.63 ±0.03 -5.12 ±0.07 



Note. — Combined fit used a single parallax parameter for both maser 
spots relative to the two background sources; a single proper motion was 
fit for each maser spot relative to all two background sources. 
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Table 3. W51 Parallax & Proper Motion Fit 



Maser Vlsr Background Parallax jji^ fj,y 

(km s~^) Source (mas) (mas y~^) (mas y~^) 



56.0 J1922+1504 0.240 ± 0.078 -2.54 ± 0.09 -5.45 ± 0.33 

56.0 J1922+1530 0.055 ±0.111 -2.33 ±0.13 -5.49 ±0.17 

56.0 J1924+1540 0.377 ±0.149 -2.60 ±0.20 -5.42 ± 0.17 

56.4 J1922+1504 0.155 ±0.127 -2.51 ± 0.15 -5.56 ±0.21 

56.4 J1922+1530 0.023 ± 0.062 -2.33 ± 0.12 -5.60 ±0.05 

56.4 J1924+1540 0.317 ±0.094 -2.59 ±0.13 -5.52 ± 0.10 

56.0 combined 0.195 ±0.071 -2.49 ± 0.07 -5.45 ±0.14 

56.4 -2.48 ±0.08 -5.56 ± 0.08 



Note. — Combined fit used a single parallax parameter for both maser 
spots relative to the three background sources; a single proper motion was 
fit for each maser spot relative to all three background sources. 



